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Abstract 

Q \ We study the possible tetraquark interpretation of light scalar meson states ao(980), /o(980), k, 

o 

CsJ ' a within the framework of the non-relativistic potential model. The wave functions of tetraquark 

■ states are obtained in a space spanned by multiple Gaussian functions. We find that the mass 
^ i spectra of the light scalar mesons can be well accommodated in the tetraquark picture if we 

' introduce a three-body quark interaction in the quark model. Using the obtained multiple Gaussian 

■ wave functions, the decay constants of tetraquarks are also calculated within the "fall apart" 
^1 mechanism. 

PACS numbers: 14.40. Cs, 12. 39.Pn,13.25.Jx 

> 

in 

d 

o 

> 

X 
J3 



Electronic address: dwz@th.phy.pku.edu.cn 



I. INTRODUCTION 



Tetraquarks were proposed decades ago. Early in 1977, Jaffe make a calculation using the 
color-magnetic interaction in the bag model [4]. He suggested that the light scalar mesons 
below 1 GeV, ao(980), /o(980), k and cr, be interpreted as a nonet of light tetraquarks. 

In recent years, the light scalar mesons are observed in decays of charmed mesons. The 

nn n 

a is observed as a peak in decay tt+tt 7r+ [3|, |J] and k in tt+tt tt^ [5]. From 

the process J/\l' — tuvr+vr", BES collaboration determined the pole position of a to be 
M - iT/2 = (541 ± 39) - i{252 ± 42) Mevje]. Also BES collaboration found a k like 
structure in the decay J/\E' ^ A'*K+7r~|3|. The accumulation of experimental data allows 
us to study the structure of the light scalar nonet based on their decay properties j^, 9, 10|. 

As a many-body system, a tetraquark state is quite different from a baryon or a conven- 
tional qq meson. The color structure is no longer trivial. It is quite sensitive to the hidden 
color structure of QCD interaction. A tetraquark state, if its existence is confirmed, may 
provide us important information about the QCD interaction that is absent from the ordi- 
nary baryons or the qq mesons. For instance, some authors had investigated the tetraquark 
system with the three-body qqq and qqq interaction, whose existence has no direct effect 
on the ordinary hadron states jlll . Il2l . Il3| . The newly updated experimental data can shed 
more light on the relation of the possible tetraquark states and QCD interaction. 

In this article, we will study the possible tetraquark state within the framework of the 
non-relativistic potential model. We will calculate mass spectra and wave functions of the 

n 

light tetraquark using the Bhaduri potential[l^. To fit the experimental masses, we will 
extend the model with the three-body qqq and qqq interaction. Using the wave functions of 
tetraquarks, we will determine the coupling constants of tetraquarks to mesons under the 
"fall apart" mechanism. 

The article is organized as follows: In Sec. [Tll we introduce the model Hamiltonian and the 
multiple Gaussian function method which is used to obtain the tetraquark wave functions. 
In Sec. mil we present the "fall apart" decay calculation with tetraquark wave function. In 
Sec. \IV\ we present the numerical results. Finally we will give a brief summary. 
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II. HAMILTONIAN AND WAVE FUNCTIONS 



In a non-relativistic quark model, usually the potentials are limited to the two-body 
interaction, which mainly consists of two parts: the A'^ ■ color interaction including the 
confinement and the Coulomb interaction of one-gluon exchange, and the A'^ ■ A^cr ■ a color- 
magnetic interaction. The Hamiltonian reads 



2m,- 



4 



(1) 



i<j 



where rrii are quark masses, Ff = -j-, and Vij is the distance between quark i and quark j. 

Among the various potential forms used in different quark models, the Bhaduri potential 
[l^ is rather simple and gives a unified description of conventional hadron spectroscopy. 

3, 0, Q, Q. The potential 



Also it is often used to discuss the tetraquark system |l3l . 
reads 



c 



\- — - D 



The parameter values are 

K = 102.67Mevfm, = 0.0326(MeV-^fm)i D = 913.5MeV, tq 
md = 337MeV, = 600MeV, = 1870MeV, nib 



0.4545fm 
5259MeV. (2) 



In a tetraquark, some new interactions which have no direct effect on the ordinary hadrons 
may have significant contribution. For instance, one can introduce the following three-body 



qqq and qqq interactions 
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V,,,ir,, r„ n) = d'^'^F^F^F^*U, exp[-(rj + + rl)/rll 
V,,,{r,, r„ nO = d^"^FtF;*F^*Uo exp[-(4 + r% + rl,)/rl]. 

In this article, since we will only treat the ground states of tetraquark, the spatial dependence 
of the three-body interaction is less important. So we will only add the following simplified 
interaction into the model Hamiltonian ([T]) 

= Uo{d'''"F;'F'F^* + d'''''F,^F'*F^*). (3) 



This interaction is diagonal in the diquark-anti-diquark color base of tetraquark 

20 



(Me Me I Vsb I [qqhlqqh) 



(4a) 
(4b) 
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An immediate consequence is that this three-body interaction has no direct contribution to 
any meson-meson couphng channel. 

To explain our calculation method, we first define some convenient coordinates for 
tetraquark system as illustrated in figure [T] 17|] , 



X3 



ri - r2, 
miVi + m2r2 

mi + 7712 



(5a) 
(5b) 
(5c) 



yi = ri - rg, (6a) 

y2 = r2-ri, (6b) 

rriiri + Triors m2r2 + m4r4 . , 

ys = (6c) 

nil + f^s ^2 + ^4 



zi = ri - r4, (7a) 

22 = r-2 - r-3, (7b) 

rriiri + m4r4 m2r-2 + m^rs , 

Z3 = (7c) 

mi + m2 + 




(a) (6) (c) 

FIG. 1: Three ways to define the relative coordinates for a tetraquark system. Filled and open 
circles represent quarks and anti-quarks respectively. 

The base wave function for the tetraquark will be the product of color, spin, fiavor and 
spatial wave functions. The color and spin SUc{3) ® SUs{2) base function we use is of the 
following diquark anti-diquark coupling form: 
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• s = 







ai = |3l2334)c ® |0l2034)^ 
"3 = |6l2634)c® IO12O34). 



"2 = |3l2334)c ® |ll2l34); 
"4 = 1612634)0 ® 1 112134); 



(8) 



S 



ai = |3i2334)c (S) |0i2l34)s, 
"3 = |3l2334)c ® |ll2l34)s, 
"5 = |6l2634)c ® |ll2034)s, 



"2 = |3l2334)c ® |ll2034); 
04 = 1612634)0 ® |0l2l34); 
"6 = 1612634)0 ® |ll2l34); 



(9) 



S = 2 



"1 = |3l2334)c ® |ll2l34); 



"2 = 1612634)0 ® |ll2l34); 



(10) 



Here the color wave function of the two (anti-)quarks is labeled by SUc{3) dimension and 
the spin wave function is by the total spin. 

The anti-symmetric diquarks [ud], [us], [ds] form the 3 representation of flavor SUf{3). 
The 3 diquarks and 3 anti-di quarks further form a tetraquark nonet. They are assumed to 



be the light scalar mesons 



19j . So the flavor wave functions are: 



ao(/ = 1, J3 = 0) = -^{[us][us] - [ds][ds]) 

f,(l = 0) = l=([us][us] + [ds][ds]) 

ao{I = 0) = [ud][ud] 
k'^ = [ud][sd] 



(11a) 

(lib) 

(11c) 
(lid) 



As for the spatial wave functions, we will start from the multi-dimensional Gaussian 
function 

g'{xi, X2, xs) = exp - ^ A'-jXi ■ Xj ] , (12) 



where are the function parameters. The wave function of this form is well convergent 



and there exists many analytical expressions 



"or different matrix elements. We will use it to 



construct the spatial wave function base 
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2l| 
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Under the hypothesis of Jaffe, the color-spin wave function of "good" diquark is the 
symmetric one |3i2)c ® |0i2)s. As the flavor wave function of the scalar nonet state is anti- 
symmetric, so the spatial wave function should be symmetric. That is the spatial wave 
function of tetraquark state should be invariant under Xi —>■ —Xi and/or X2 —X2. If we 
use the Gaussian function (IT^ as the base wave function, it is easy to see that [17j 

A12 = A23 = A31 = 0. 

We will use the following symmetric combination as the base function 

^p'ixi, X2, X3) = ^[g'{xi, X2, X3) + g'{-xi, X2, x^) 

+ g'i^l, -X2, X3) + /(Xi, X2, -X3)]. (13) 

If the non-diagonal parameters Aij{i 7^ j) are small, we have 

ij'{xi, X2, X3) ^ exp [-{Al^xj + A^^xl + Al^xfj] 

X [1 + 2A'^{xi ■ X2f + 2A'^{xi ■ xsf + 2Af^[x2 ■ x^f\ . (14) 

This allow us to study the correlations in the quark alignment. 

We will choose n independent symmetric Gaussian functions (|T3l) . s = l,2,...,n, to 
span an n-dimensional nonorthogonal basis. The n independent Gaussian functions are 
obtained by the following process. First, we use one such symmetric Gaussian function as 
the test wave function in variation to determine a base parameter set Aij. The matrix (Aij) 
will be specified by three principal values denoted ' ^422^ ^433^ and three Euler angels 
(a,/3,7) which specified the orientation. Then a complete parameter set A^As = 1,2, ...,n) 



is generated by first scaling to the principal values [17] 

Af^ = (15) 

where = —k, —k + 1, k — l,k, {2k + 1)^ = n, and d a scaling factor; Then we make an 
Euler rotation (a,/3, 7). 

By diagonalizing the Hamiltonian in the above nonorthogonal basis, we will obtain the 
mass and wave function of tetraquark states. The wave function can be expressed in the 
above base functions as 

\T)=(j)fJ2Cisair (16) 



where 0/ is the flavor wave function, Cis is the superposition coeflicient. 

Similar to case in pseudo-scalar mesons, the 1 = members /o, ctq in the scalar nonet 
will mix with each other. To consider the mixing, we further introduce a mixing angle (j) jo] 

/ = /o cos (/) + (Jo sin 0, (J = — sin(/)/o + cos</)cro- (17) 

Then / and a are the physically observable states. In this article, we do not discuss the 
underlying mechanism of this mixing. So we will merely treat the mixing angle as one 
additional parameter. 



III. DECAY PROPERTY OF TETRAQUARK STATE 



Several authors have used the effective Lagrangian with SUf{3) symmetry to discuss the 



decay of light scalar nonet 



Here we can calculate the coupling constants using the 
tetraquark wave functions. The general coupling Lagrangian reads 

TT ■ TT 



C = fo 



gfoKKKK + ... 



2 

vr ■ vr - 



+ 9KKn {KtK ■ TT + h.c) + ... 



(18) 



At present, the quark interaction underlying those meson decaying couplings is still un- 
clear to us. Here we will assume that the decaying is the fusion process and can be depicted 
by the "fall apart" mechanism in flgure O More speciflc, we assume that the coupling con- 
stant of a tetraquark T to two mesons Mi and M2 is proportional to the wave function 
overlapping 

gTMM oc (M1M2 I T). (19) 

The meson wave functions will also be approximated by multiple Gaussian wave functions 
determined by a similar variation process 



\M)r = <PsY,Csr{r), 

s 

where 0/ is meson flavor wave function, and the spatial base function is 

^,(r) = 



(20) 



(21) 



A tetraquark system gig2Q'3Q'4 can fall apart into two different flavor combinations qiq^ + q2<lA 
and qiqi + ^2^3, and the corresponding final meson-meson states are different 

\MiM2)i = \Mi)y,\M2)y,, (22) 
|MiM2)2 = |Mi),jM2).,. (23) 

The spatial wave functions are in the coordinates yi and Zi defined in Eqs. ([6]) and ([7]) 
respectively. 



T . 




FIG. 2: "fall apart" mechanism for decays of q q tetraquark states 

In the decay of the light scalars to pseudo-scalar mesons, we need to consider the rj-r]' 
mixing 

T] = COS Orjq + sin 9r]s, (24a) 
rj' = — sin 9riq + cos Orjs, (24b) 



where r]g = -^{uu+dd), rjg = ss and sin 6 = —0.608 [22|. We obtain the following expressions 
for the coupling constants (a proportional constant is dropped) 

9 h->m = sin cos OA /o-.r,,+r,, (25a) 

9fo^m' = ~ sin^^)^/o-7y,+r,, (25b) 

9fo^v'v' = - sin cos 0Af^^.^^+.^^ (25c) 

9fo-^KK = -^^fQ^K++K- (25d) 



8 



9(TQ—*TT-K 2 ^O"0— >7r + +7r~ (26a) 

9ao^m = \ cose'^A^„^r,,+r,, (26b) 

gao^m' " ~ ^ ^^<xo^'?,+'?, (26c) 

^^0^^/^, = i sin ^^,,0^^,+^, (26d) 

9a^^v = ^ sin (27a) 

^a--r,' = ^ COS dA^+^^+^^^ (27b) 

ga->KK = -^Aa+^K++KO (27c) 

Qk-^ttK = — v4^+^^++xo (28a) 

5f«,^^if = ^ cos 9A^+^^^ (28b) 

= sin 6'A«,+_,;^++^^ (28c) 



Besides the explicit flavor overlapping factors, At^mm is the overlapping of the color, spin 
and spatial wave function. 

After considering the a-fo mixing effect. The coupling constants qt^mm for the decays 
of a and /o are further modified to 

9f^MM = COS (j)gf„-,MM + sin (f)gao-*MM, (29a) 

ga~*MM = - sin (pQfo-^MM + COS (pQao-^MM (29b) 

IV. NUMERICAL RESULTS 

Bhaduri potential gives a unified description of the spectroscopy of ordinary mesons and 
baryons. The Hamiltonian ([T]) itself is an eigenvalue problem of the differential equation 
which can be solved numerically. However, the multiple Gaussian function method can still 
give an impressively good approximation of the ground state mesons and the Gaussian wave 

9 



function is rather simple to use. In Table HI we show some results of the pseudo-scalar meson 
calculation. We can see that the multiple Gaussian function method greatly improve the 
single Gaussian function approximation. 

= m^, (MeV) 136 250 137 

rriK (MeV) 520 582 521 

m^^ (MeV) 758 800 758 

TABLE I: Pseudo-scalar meson calculation. In col. 1, we direct solve the Schrodinger equation. In 
col. 2, we use the variation method with a single Gaussian function. In col. 3, we use the multiple 
Gaussian function method with 7 Gaussian functions. 

Now we turn to the tetraquark calculation. In our calculation, the scaling factor is fixed 
to be d = 2. We will take = 1, i.e., the wave function space is spanned by 3^ = 27 
Gaussian functions. In the light scalar tetraquark, as we assume that the flavors of diquark 
and anti-diquark are antisymmetric and the spatial wave function is symmetric, so the color 
and spin wave function must be the symmetric ai and in eq. ([H]). 

First, we will consider the original Bhaduri potential without the three-body quark in- 
teraction ([3]). We obtain the following masse values 

M„ = 687MeV, M« = 1067MeV Ma, = Mj^ = 1371MeV. (30) 

We can see that the masse values are about 300 MeV higher than the experimental values. 
We can calculate the possibility of a tetraquark state in different color-spin structure 

r ^ 

= l[dx,\{a,\^)f . (31) 

^ i=i 

The color-spin contents of the tetraquark nonet without three quark interaction are presented 
in Table [TTl We see that the color content is mainly the 6x6 component which disagrees 
with Jaffe's "good" diquark hypothesis. 

Next, we will turn on the three body interaction with Uq = 0.333GeV. We find that light 
tetraquark masses are 

M„ = 443MeV, M« = 744MeV M„„ = Mf^ = 985MeV, (32) 
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p 



a 
0.30 
0.70 



K 

0.30 
0.70 



ao> /o 
0.29 
0.71 



TABLE II: Contents of tetraquarks without three-body interaction. 



which are in agreement with the experimental values 23l |: 



M„ = 800 ± 400MeV, 
Ma, = 984.7 ± 1.2MeV, 



= 840 ± 80MeV, 
Mf, = 980 ± lOMeV. 



(33) 



The color-spin contents o: 
"good" diquark picture j^, 



the nonet are shown in in Table IIIII which also agree with the 



24|. 



ai 



Pr 



a 
0.80 
0.20 



K 

0.88 
0.12 



ao> /o 
0.92 

0.08 



TABLE III: Contents of tetraquarks with three-body interaction. 

In our calculation, the tetraquark wave function is symmetric under the coordinates 
reflections Xi — > —xi and/or X2 —X2. It is easy to see that the expectation values 



(34) 



a/ {x\) and ^J {x"^ are the radii of mean square (RMS) of diquark and anti-diquark respec- 
tively. The quark and anti-quark RMS in tetraquark is 

.2 _ (Eto"^^(^^ - RcmY) 



(Jil2{xl) + ^lZi{xl) + Hi2,m{xI) 

rui + m2 + ms + rrii 



(35) 
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where 



-KCM — -=^4 , 1-30 j 

^^ij = (37) 

rrii + mj 

_ (mj + mj)(mfc + m«) 

- ^ : : • [■io 

rui + TUj +mk + mi 



The RMS values are tabulated in Table fllVD. 



a K ao, fo 

^/(xf) 0.70 0.72 0.70 

^/(x^ 0.70 0.69 0.70 

^/{x^ 0.54 0.58 0.56 



TABLE IV: The RMS values in fm. 

However, the spatial wave function (fT3|) is beyond the usual tetraquark assumption. 
Usually a tetraquark is assumed to be constructed from the "good diquark". The inner 
orbital angular momentum of the (anti-) diquark in a tetraquark is zero. So the relative 
angular momentum between the diquark and anti-diquark in the scalar tetraquark is also 
zero. The spatial wave function will has the form 

i/j{xi,X2,X3) =ip{xl,xl,xl), (39) 

i.e., all the Xi are in S-waves. Our choice f|T3l) is beyond the above assumption (this can be 
easily see from eq. (fT4l)). If eq. (!39ll holds, we will have following identity 

{{x.-x,r) = l{x^){x'^) (z^j). (40) 



We will measure the deviation from (1391) of a tetraquark state by 

_ 3{{xi-Xjy) 

The numerical ejj values are listed in Table El The small nonzero e values means that the 
tetraquark states are indeed not pure in S-wave. There is always some D-wave mixing. 

With the obtained wave functions, we can calculate the wave function overlapping 
in Eqs. ( I25|l - (!28l) to get the coupling constants. The results are collected in Table |Vll 
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a 


K 




ao> /o 


ei2 


0.14 


0.23 




0.08 


ei3 


0.21 


0.22 




0.13 


£23 


0.21 


0.24 




0.13 




TABLE V: The eij 


values of tetraquark wave function. 






(To — > vr"*" + TT^ 


K+ ^ 7r+ + 


^ K+ + dd a+ ^ K+ + 


a+ ■ 


7r+ + 


10.75 


9.37 


9.37 8.16 




8.38 



TABLE VL Tetraquark-meson-meson wave function overlapping in color, spin, spatial space 
At^m^M2 (unit GeV~3/2). 

According to ref. jo!, the scalar isoscalar mixing angle (p in eq. f[T71) will be fixed by the ratio 
d'f^Kx/ d'f^TrTT ~ ^■'^^ with eq. fl29|) . This gives = 16.8°. The ratios of coupling constants 
for scalar meson decays are listed in Table Ivnl Similar to Bugg's calculation |9|], although 
most of the experimental ratios can be fitted within a factor 2, g'j^^^/ g'j_^^^ is far above 
experimental value. 



V. SUMMARY 



In summary, we have performed a tetraquark calculation of light scalar mesons using 
the quark potential model. If we only consider the two-body quark interaction as in the 
conventional hadron calculation, the masses of the tetraquark states will be several hundred 
MeV higher than experimental data. Also the major component of the light tetraquark wave 
functions consists of the color sextet diquark and anti-diquark. After including a three-body 
interaction in the Hamiltonian, we find the masses of the light tetraquark nonet agree with 
experimental data and their wave functions are composed of mainly the "good" diquark and 
anti-diquark. We have used the multiply Gaussian function to approximate the tetraquark 
wave functions and noticed that there is a small mixing of D-waves in the wave functions. 
With wave functions we obtain, we also calculate the coupling constants for scalar meson 
decays according to the "fall apart" mechanism. By introducing the isoscalar mixing angle 
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Analysis of ref . 2^ Analysis of ref. 9 


our results 


Exptte, 26, 27, 281 


gin 






n fin 


U.4U =b U.Uo 


0.39 


0.75 ±0.11 








A 01 


41:. Zl Zt U.4I:D 


4.21 


4.21 ±0.46 


^/-* 


KK' ^ao- 




z.zo 


n no 1 n m 

u.yo ± o.oi 


0.92 


2.15 ±0.4 






— 'TTTj 


U.io 




1.71 












i.U / it u.io 


1.15 


< 0.33 




■*KKl 


->7r7r 


4.8 X 10"'^ 


0.03 ± 0.01 


0.04 


0.6 ± 0.1 


i/cr- 




>7r7r 


0.05 


0.23 ±0.02 


9^ 


n on _|_ n 04 


5i 


-^ttkI 9I- 


->7r7r 


0.78 


0.58 


0.83 


(2.14 ±0.28) to (1.35 ±0.10) 


9l 


■*r]Kl 9\- 


>ttK 


0.12 


0.20 ±0.01 


0.21 


0.06 ± 0.02 


9l- 




■*nK 


0.006 


0.13 ±0.01 


0.12 


0.29 ±0.29 



TABLE VII: Ratios of coupling constants for light scalar meson decays, with (p = 16.8°. 

0, we obtain a fit of the ratios of coupling constants for scalar meson decays similar to other 
analysis based on tetraquark picture. 
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